Mixture of bosonic and spin-polarized fermionic atoms in an optical lattice 
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We investigate the properties of trapped Bose-Fermi mixtures for experimentally relevant param- 
eters in one dimension. The effect of the attractive Bose-Fermi interaction onto the bosons is to 
deepen the parabolic trapping potential, and to reduce the bosonic repulsion in higher order, leading 
to an increase in bosonic coherence. The opposite effect was observed in 87 Rb - 40 K experiments, 
most likely due to a sharp rise in temperature. We also discuss low-temperature features, such 
as a bosonic Mott insulator transition driven by the fermion concentration, and the formation of 
composite particles such as polarons and molecules. 
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I. INTRODUCTION 

Systems of interacting bosons and fcrmions occur fre- 
quently in nature. Usually, the bosons act as carri- 
ers of force between the fermionic particles. In high 
energy physics, quarks exchange gluons via the strong 
force, while in solid state physics electrons can interact 
via light or lattice vibrations. Most prominent examples 
of such systems are conventional BCS superconductivity 
(caused by an effective attractive interaction between the 
fcrmions induced by the electron-phonon coupling), the 
Peierls instability (a charge density wave) and the for- 
mation of polarons, which in solids are electrons dressed 
by a cloud of phonons. There are only a few condensed 
matter systems in which the influence of fermions onto 
bosons has been investigated. One of them are mixtures 
of bosonic 4 He and fermionic 3 Hc, in which a shift of the 
transition temperature between normal and supcrfluid 
4 He as a function of 3 He concentration was observed. 

In the field of ultracold gases fermions and bosons 
are on an equal footing. The choice of different atomic 
species [l|, 0,0, HI , the use of Feshbach resonances @, H, 0] 
and optical lattice potentials [1, Q give almost unre- 
stricted access to all parameters of these systems, offer- 
ing the possibility to study the influence of the species 
onto each other and to investigate open questions from 
other areas of physics in a new context. Theoretical ap- 
proaches 0, El, Ei El EI EI M EE El, El have 

proposed a whole variety of quantum phases present in 
homogeneous Bose-Fermi mixtures at low temperature, 
ranging from a charge-density wave, over a fermionic 
pairing phase, to polaronic properties, and even to phase 
separation. 

In recent experiments two groups independently suc- 
ceeded in the stabilization of bosonic 87 Rb and fermionic 
40 K in a three-dimensional optical lattice [1,0. They fo- 
cused on the loss of bosonic phase coherence and on the 
increase of the bosonic density by varying the fermionic 
concentration. The trapping potential and the finite tem- 
perature make the interpretation of the observed quanti- 
ties however challenging. 



Here, we study the interplay of a trap, finite temper- 
ature and strong intcrparticlc interactions, which lead 
to physics quite different from the homogeneous case. In 
particular, while the addition of fermions induces a quan- 
tum phase transition from the Mott insulating to the su- 
perfluid phase at larger bosonic repulsion strength in a 
homogeneous lattice, the presence of a trapping potential 
makes the situation more involved because of an extra, 
effective strongly-inhomogencous trapping potential. De- 
spite the large Bosc-Fcrmi coupling, it turns out that our 
results for the trapped, mixed system can be well under- 
stood in terms of first and second order corrections to 
the bosonic Hamiltonian. For low bosonic and fermionic 
densities, we illustrate that the formation of bound pairs 
invalidates the picture of a perturbational correction by 
the fcrmions on the bosons. 



II. MODEL 

A mixture of bosonic and spin-polarized fermionic 
atoms in an optical lattice can be described by the Bose- 
Fcrmi Hubbard Hamiltonian, 
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where cf(&f) and Ci(bi) arc the corresponding creation 
and annihilation operators for the fermions (bosons), and 
nx,i is the number operator on site i for species X=B,F. 
The Jx-terms are the hopping, the ex-terms describe the 
external trapping potential, and Ubb and Ubf denote 
the on-site interaction strength between bosonic atoms 
and bewcen a fermionic and bosonic atom, respectively. 
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FIG. 1: (Color online) (a) Shift in critical U BB of the Mott 
transition for attractive Ubf- The system consists of 13 
fermions and 64 bosons on a homogeneous lattice of 64 sites. 
Both species have unit hopping and the inverse temperature 
is ft = 64/ J B . The critical value U B b/Jb = 3.28 ± 0.04 at 
Ubf/Jb = is taken from Ref. [2^]. At finite Ubf the tran- 
sition is located where the Green function has the same al- 
gebraic decay as in the purely bosonic case. The solid curve 
is a parabolic fit. (b) Comparison between the calculation 
in the presence of fermions and their approximation by a 
site-dependent potential for a trapped system of 60 sites, 40 
bosons, 8 fermions, j3 = 1/Jb and optical potentials Vo = 6-Er 
(Ubb/Jb = 11.89). Here, Eb. = h 2 k 2 /2m,B,b is the bosonic re- 
coil energy. 71b(f) denotes the bosonic (fermionic) density for 
the mixture, no is the density obtained in the approximation. 
Analogous for the density fluctuations Kb and kq. 



The effective parameters of the Bose-Fermi Hubbard 
model are deduced from the experimental parameters of 
Ref. [1, [j| using a tight-binding approximation [ill, fpoll - 
Taking the scattering len gths as aeB/ao = 102 ± 6 [2ll | 
and clbf/o-o = —205 ± 5 [22j, where ao is the Bohr ra- 
dius, we note that Ubf/Ubb ~ —2, a ratio which is al- 
most constant for all optical lattice depths. We took a 
wavelength A = 1064nm for the optical lattice potential, 
and frequencies lob = 27r • 30Hz and cup = 2n ■ 37Hz for 
the harmonic confinement. To determine the state of 
the mixture we use two numerically exact methods: at 
finite temperatures the canonical two-body Bose-Fermi 
worm Quantum Monte Carlo (QMC) algorithm [13], and 
at zero temperature the density-matrix renormalization- 
group method (DMRG) 



III. EFFECTS OF FERMIONS ON BOSONS 
A. Induced potentials and interactions 

In a mixture, the lowest order effect of the fermions is 
a mean-field shift [/RF(iF.i) of the potential experienced 
by the bosons [HI, Hj] . In a homogeneous system with 
a fixed particle number the shift in the potential has no 




FIG. 2: (Color online) Dependence of bosonic (red, upper 
curve) and fermionic density (blue, lower curve) profiles on 
the interspecies interaction strength. In the system there are 
60 sites, 50 bosons, 20 fermions, optical potential is Vq = 3Eb. 
(Ubb/Jb = 4.26), and the inverse temperature (3 = 4.26/Jb- 



consequences besides adding a constant to the energy. 
The next order effect is an induced attractive interaction 
between the bosons [H|,[24|] ■ A similar effect is well known 
from conventional superconductivity, where the phonons 
(bosons) induce an effective electron-electron interaction. 
The induced interaction shows up most clearly in a shift 
of the critical Uq B /Jb of the bosonic superfluid-Mott 
transition while varying the interspecies interaction as 
shown in Fig. []Ja) . As expected for an induced attractive 
interaction, we find a shift to larger values of Ubb/Jb- 
At small Ubf the shift is proportional to U BF in agree- 
ment with perturbative calculations (24|. Therefore the 
presence of fermions can induce a phase transition from 
a Mott-insulating to a supcrfluid phase. 

In the following we discuss how a parabolic trap and 
finite temperature change this picture. To separate the 
effect of the effective trapping from the induced interac- 
tion, we generate an effective site-dependent potential for 
the bosons by replacing the Bose-Fermi interaction op- 
erator by the effective potential /LtjfiB,j = E^bf (^f, i)^B,z ■ 
This deviates from the disordered chemical potential ap- 
proach of Ref. Q and from a mean-field approximation: 
The exact fermionic density distribution of the mixture 
serves as input for a second, purely bosonic simulation. 
In Fig. [Hh), we compare the resulting bosonic density 
and compressibility profiles. We observe that the den- 
sity profiles are quite well reproduced, confirming that 
the dominant effect of the fermions is to modify the ef- 
fective potential for the bosons. However, looking at 
higher order quantities such as the density fluctuations 
( n B i) ~ ( n B,i) 2 we find significant discrepancies. In par- 
ticular, we see Mott plateaus in the approximation (sig- 
naled by dips in the variance of the density in Fig. [lib)) 
that are absent in the full QMC simulation. Around 
these dips the difference between the two curves is around 
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eighty percent. This is a clear signature for a fermion- 
induced attractive Bose-Bose interaction, reducing the 
bare repulsion Ubb (cf. Fig.QJa)). We note that the vis- 
ibilities, discussed below, are rather well reproduced in 
the approximation, indicating that the effective potential 
is the dominant effect as far as this experimental quantity 
is concerned (which is surprising seen the large values of 
Ubf)- 

Having gained an understanding of the relevant mech- 
anisms, we proceed in section llll Bl with the results of two 
simulations where we vary experimental control parame- 
ters, namely the Bose-Fcrmi coupling and the fermionic 
concentration, and then compare our results to experi- 
ment in section llll CI 
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B. Results at low and constant temperature 

Fig. [2] shows density profiles for different values of the 
attractive interaction strength Ubf- In the absence of a 
boson-fermion interaction, all particles are smeared out 
over the lattice. Turning the interspecies interaction on, 
we see in Fig. [2] that both species accumulate in the trap 
center. The fermions are pinned down in the trap center 
(cf. Fig. [2]), despite their light bare mass. They lower the 
effective potential in the center of the trap as C/BF^F.i)) 
causing the accumulation of bosons. 

The effect of an inhomogencous effective potential can 
be seen even more clearly by varying the fermionic con- 
centration instead of the interaction strength Ubf- In 
Fig. [3] we show the dependence of the bosonic visibility 
on the number of fermions with a fixed number of bosons, 
a setup similar to recent experiments [H, Q . The bosonic 
visibility is defined by v = (pb(0) — Pb(.tt))/{pb{0) + 
PbM) where Pb(^) is the value of the bosonic momen- 
tum distribution at momentum k. The visibility is often 
taken as a measure of the coherence of the bosons. 

For shallow lattice potentials, a slight increase in the 
visibility for an intermediate number of fermions is seen, 
due to an increase in the bosonic density in the center of 
the trap caused by the fermions. For moderate lattice po- 
tentials one observes a complex, non-monotonic behavior 
with large variations. If a few fermions are admixed to 
a bosonic system, the fermions - spreading over several 
sites in the center of the trap - cause a strongly inhomo- 
geneous effective potential for the bosons. The effective 
potential exhibits a deep minimum in the center of the 
trap and causes the bosons to accumulate in this region. 
If the purely bosonic system was superffuid (cf. Fig. [3j 
Ny- = 0, Vo = 3 .Eft) the effective potential causes a su- 
pcrfluid state with higher filling in the center of the trap, 
slightly increasing the bosonic visibility. If the bosonic 
system exhibited a broad jib ~ 1 Mott plateau (cf. inset 
Fig. [3J Np = 0, Vo > 6£ft), this plateaux is partially 
destroyed resulting in a large rise of the visibility. The 
mechanism holds until there arc enough fermions present 
to form a band insulating region (cf. Fig. [3J Np « 14). 
For such and higher fcrmion numbers, the effective potcn- 



FIG. 3: (Color online) The bosonic visibility v as a function 
of fermionic number for a system of L = 60 sites, TVs = 40 
bosons and inverse temperature /3 = 1/Jb for optical po- 
tentials Vo = 3,6,8£ R (or Ubb/Jb = 4.26,11.89,21.58, re- 
spectively). The inset shows the bosonic density profiles at 
Vo = 6-Er for N f = 0, 40, 4, 20, 16 bottom to top in the trap 
center. 




FIG. 4: (Color online) Change in bosonic visibility when in- 
verse temperature (3 is increased for a system of 100 sites, 50 
bosons and 30 fermions. 



tial induced by the fermions follows the curvature of the 
external trap over the region occupied by the fermions 
with a sudden increase at its boundaries. The number 
of fermions sets the length of an effective system for the 
bosons, and controls the bosonic filling. In this approxi- 
mately parabolically trapped effective systems insulating 
regions can form for Nb/Np = 40/20 (as shown), but also 
for N B /N F fn 60/30 or 20/10, yielding strong dips in 
the visibility [2(1 12^. 
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FIG. 5: (Color online) (a) Fermionic (pf(A;)) and molecu- 
lar Bose-Fermi (pp{k)) momentum profiles for the same sys- 
tem. Error bars are shown, but are smaller than the point 
size, (b) Bosonic and fermionic ground state density pro- 
files for a system with parameters L = 60,iVB = Nf = 8, Vo = 
4-Er(C/bb/ Jb = 6.13) determined using DMRG (QMC results 
agree within error bars). The bosonic and fermionic density 
profiles almost coincide. 



C. Comparison with experiment 

In apparent contradiction to our low-temperature pre- 
dictions of an increase in the bosonic visibility for most 
numbers of admixed fermions (in Fig. [3|), the experi- 
ments [H, [3] show a strong decrease for moderate values of 
the lattice potential. Assuming that entropy is conserved 
when ramping up the lattice, the effective temperature of 
a 87 Rb - 40 K mixture rises dramatically because of the dif- 
ferent temperature dependence of fermionic (Sf oc T/Tf 
for an ideal gas) and bosonic (Sb oc (T/T c ) 3 for an ideal 
gas) contribution to the entropy Q. In Fig. [4] we increase 
the temperature (decrease (3 — l/fc^T) at fixed opti- 
cal potential and atom number and find a drastic drop 
in the visibility as the temperature is increased around 
(3 w 2/Jb- A rise in the temperature of the bosonic 
atoms in the presence of fermions can cause a large de- 
crease of the bosonic visibility and is thus the most likely 
explanation for the experimental results. 

IV. AT LOW DENSITIES 

We finally discuss the physics at low densities. Bound 
pairs ("molecules") of one boson and one fermion can 



now be formed for moderately deep optical lattices. 
First signs of this pairing can be seen in the two- 
body Bose-Fermi momentum distribution shown in Fig. O 
(a). The momentum distribution of these fermionic 
molecules shows a sharper Fermi edge compared to the 
bare fermion. It also tends to zero at larger momenta, 
showing that these molecules are a better description of 
the system than the bare fermions and bosons. The for- 
mation of molecules is also well supported by coinciding 
charge modulations in the bosonic and fermionic densi- 
ties (Fig. Efb)). For the parameters chosen in Fig. [5fb) 
the density modulations are Friedel oscillations due to 
trap, but for larger lattice depths a density wave can be 
formed [13, El- 



V. CONCLUSION 



In conclusion, we have simulated the trapped one- 
dimcnsional Bose-Fermi Hubbard model. The interplay 
between temperature, trap, optical potential and parti- 
cle number is very rich and non-universal. The domi- 
nant effect is the creation of a strongly inhomogeneous 
trapping potential by the fermions. In higher order the 
fermions induce an attractive interaction between the 
bosons, which should lead to an increase in the bosonic 
visibility. However assuming a rise in temperature when 
ramping up the lattice, a decrease in the visibility is 
found analogous to the experimental observation in the 
87 Rb - 40 K samples. If temperature remains low, one 
could observe a Mott transition driven by the fermionic 
concentration, and observe the formation of molecules. 
The same effects are expected for higher dimensions, 
since the underlying mechanisms do not depend on di- 
mensionality. Our study explicitly shows that the effects 
of temperature, particle number, adiabatic processes and 
trapping potential have to be taken into account carefully 
when analyzing cold-atom experiments. 
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